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Abstract 

Counterstreaming plasma structures are widely present in laboratory experiments and astrophys- 
ical systems, and they are investigated either to prevent unstable modes arising in beam-plasma 
experiments or to prove the existence of large scale magnetic fields in astrophysical objects. Fila- 
mentation instability arises in a counterstreaming plasma and is responsible for the magnetization 
of the plasma. Filamentationally unstable mode is described by assuming that each of the counter- 
streaming plasmas has an isotropic Lorentzian (kappa) distribution. In this case, the filamentation 
instability growth rate can reach a maximum value markedly larger than that for a a plasma with a 
Maxwellian distribution function. This behaviour is opposite to what was observed for the Weibel 
instability growth rate in a bi-kappa plasma, which is always smaller than that obtained for a 
bi-Maxwellian plasma. The approach is further generalized for a counterstreaming plasma with 
a bi-kappa temperature anisotropy. In this case, the filamentation instability growth rate is en- 
hanced by the Weibel effect when the plasma is hotter in the streaming direction, and the growth 
rate becomes even larger. These effects improve significantly the efficiency of the magnetic field 
generation, and provide further support for the potential role of the Weibel-type instabilities in 
the fast magnetization scenarios. 
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The electromagnetic instabilities of Weibel-type are driven by the velocity anisotropy of 
the plasma particles, whether it resides in a plasma with an electron temperature anisotropy 

Weibel instability [l|, Q] ) , or in a counterstreaming plasma system (filamentation instability 

3]). Davidson 0], has combined for the first time these two sources of free energy showing 
that the Weibel-type instabilities release the excess of free energy stored in "the relative 
directed motion between plasma components and/or an anisotropy in plasma kinetic energy" 

4J. Moreover, in order to describe a counterstreaming plasma with temperature anisotropies, 
Davidson [4] defined a new counterstreaming distribution function including a bi-Maxwellian 
electron temperature distribution function, and derived a dispersion relation, which admits 
purely growing solutions propagating perpendicular to the electron streaming direction. In 
a counterstreaming plasma with electron temperature anisotropies, the both free energy 
sources are present. Their contributions cumulate leading either to enhancing or quenching 
the electromagnetic instability, depending on the plasma kinetic energy perpendicular to 

;he_propagation direction of the instability, if it exceeds or not the parallel kinetic energy 

,la tz|. 



The stability properties of the counterstreaming plasma have been extensively investi- 



gated 



q, y, q y, (2 



10, 
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, and the general observation is that the effect of a finite 



plasma temperature has been revealed only for a Maxwellian electron distribution function. 
Further investigations have been carried out in order to examine these instabilities in a 
magnetized counterstreaming plasma with a bi-Maxwellian electron temperature anisotropy 
14 I la . llo l. As expected, the presence of a dc magnetic field in the streaming direction 

nnn 

inhibits the Weibel-type instabilities [8|, Hjj, uM- We therefore ignore any influence of an 
ambient magnetic field in the present investigation. 

Space plasmas are observed to possess a non-Maxwellian particle distribution function 
with a high energy tail 19|, [20|. The most common distribution function used to model such 
natural plasmas is the generalized kappa distribution [21] 
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where k is the spectral index (for typical space plasmas, k generally lies in the range 2- 
6), and thermal velocity 9 is related to the particle temperature by 9 2 = [(2k — 3)/k] v\ 
and v\ = (kBT)/m. The Maxwellian and kappa distributions differ substantially in the 
high energy tail i\v/9\ ^> 1), but in the limit k — > +00, f K degenerates into a Maxwellian 



distribution function. The above /t-distribution has been used to explain different waves 
and instabilities and their related effects in isotropic plasmas where the velocity distribution 
function is assumed to deviate from the Maxwellian thermal equilibrium due to the presence 



of high energy tails 
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Very often plasma particles have a preferred direction of motion in space (e.g. a high 
temperature direction, a stationary magnetic field direction, or a streaming direction) and 
their velocity distribution function is anisotropic. In order to describe the anisotropic plas- 
mas with high-energy tails, one can use a two-temperature bi-kappa distribution function 
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which approaches the bi-Maxwellian distribution function as k —>■ +oo. The temperature 
anisotropy is defined as A = (T y /T) — 1, with 

V Tj Oy = — V Ty , (3) 



and 



ksT _ 2 _ ksT y 



v\ = '^L. (4) 
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The bi-kappa distribution function has been earlier attributed to magnetized plasmas 



providing a useful tool for a more realistic characterization o 
celeration mechanisms in the planetary environment 



' the turbulent fields and ac- 
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261 ] . Moreover, Hellberg and 



Mace (2002) have generalized the plasma dispersion function for a kappa-Maxwellian distri- 
bution, and investigated the electrostatic waves arising in such an anisotropic plasma. The 
kappa-Maxwellian distribution function combines a one-dimensional k distribution along a 
preferred direction in space, and a two-dimensional Maxwellian in the plane perpendicu- 
lar to this direction. This choice is justified physically by the fact that in many instances 
there is equilibration and isotropization in the perpendicular plane, but there is preferential 
acceleration, leading to a power-law distribution along another preferred direction 

Recently, a bi-kappa distribution function of form (T5]) has been used to evaluate the 
Weibel instability rates in an unmagnetized plasma [2;]. For small values of k, the growth 
rate is reduced, but as n is increased, the growth rate enhances approaching the results 
obtained with a bi-Maxwellian distribution. However, it has already been shown j^, (], 
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that, in a counterstreaming plasma with bi-Maxwellian temperature anisotropies, the 
(electro)magnetic instabilities of the Weibel-type can be markedly faster when the plasma 
is hotter along the streaming direction, or they can be stabilized by a higher transverse 
plasma temperature. Moreover, there are many space plasma systems in which colliding 
plasma shells can be described as counterstreams (e.g. relative motion of fully ionized 
gaseous matters in intergalactic medium, interaction of internal or external fireballs with 
surrounding interstellar medium in gamma-ray burst sources, flares in the solar winds). 
We are therefore entitled to extend here the investigation on the aperiodic electromagnetic 
instabilities in counterstreaming plasmas with non-Maxwellian anisotropies. 
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FIG. 1: Two counterstreaming plasmas and the electromagnetic fllamentation mode with fields E 
and B, and the wavevector k perpendicular to the direction of the streams. 



In Figure [T] we fix the orientation for the counterstreams and the fllamentation mode 
we are looking for: the electric field, E || Vx,a 5 is along the streaming direction, and the 
wave vector, k _L v% 2, is perpendicular to the streams. This electromagnetic mode will be 
solution of the dispersion relation 
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where to and k are respectively, the frequency and the wave number of the plasma modes, c 
is the speed of light in vacuum, and e yy is the dielectric tensor which describes the oscillatory 
properties of the plasma system. From the linearized Vlasov-Maxwell equations, one simply 
obtains 
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where the unperturbed velocity distribution function / a o( v ) for the particles of species a 
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is considered normalized by J dv / a o(v) = 1, and 
frequency. 

We now introduce the distribution function 
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(47m a e 2 /m a ) 1//2 is the plasma 
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which describes the electron plasma counterstreams with bi-Kappa temperature anisotropies. 
This distribution function has to be added to the list of Lorentzian-type particle distribu- 
tion functions provided in Table I of Ref. [2l| . For the sake of simplicity, we assume 
symmetric counterstreams with equal densities, w p , e ,i = w p,e,2 = u P;e , and equal temper- 
ature anisotropies, A\ = A2 = A = (T y /T) — 1 , and moving with the same velocities, 
l^i I = V2 = Vq. Moreover, we assume that imobile ions form the neutralizing background, 
and the electron plasma counterstreams are considered homogeneous, and charge and current 
neutralized. 

Inserting ((7j) into ([6]) and the resultant expression into ([5]) we derive the dispersion relation 

k 2 c 2 



yy uj 2 e 2 







;i ~2^ ) + SMS) 



(8) 



where 



1 + 



p,e 



A + (A + l)-Z» 



(9) 



ke K \k9Ji 

has been obtained in Ref. [2| for an anisotropic plasma with a bi-kappa distribution function 
without vq, where the simple Weibel instability arises. Otherwise, for A = Eqs. (jHJ)-(|SJ) 
simplify to 
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which describes the filamentation instability arising from counterstreaming plasmas with 
the same isotropic kappa distribution function. The dispersion equation (jSJ) contains the 
modified plasma dispersion function, defined in Ref. [211 ]. 
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and that differs only slightly from the new modified plasma dispersion function 

Z «C/) - ^jj2 r(K-i) dx x-/ ' ^ (/) > (12) 



introduced recently in Ref. [21]. It is simple to demonstrate that the kappa dispersion 
functions from above are related to each other by 



Z ! (fl=(l + ^)z.( fl + {(l-^) 



and therefore, both of them approach the plasma dispersion function |29j] in the limit of 
k — > oo. We should also remark that, in this limit, Eq. (jSJ) becomes identical to the 
dispersion relation of cumulative filamentation-Weibel mode in a counterstreaming plasma 
with a bi-Maxwellian anisotropic distribution function (see Eq. (8) in Ref. [7|). It is 
clear now that the general dispersion relation given by ([5]) and (jHJ), allows us to extend 
the investigation of purely growing filamentation instability by including the Weibel-type 
enhancing effects due to the existence of a non-Maxwellian temperature anisotropy in the 
counterstreaming plasmas. 

An analytical expression for the aperiodic solutions, 3 (a;) = uii, $t(u) = 0, of the general 
dispersion relations (jHJ) and (jHJ), can be further obtained only by using the asymptotic forms 
of the modified plasma dispersion functions provided in Refs. [2I and for very large 
or very small arguments. Here, we keep the accuracy of a general approach and evaluate 
exactly the numerical growth rates for several representative cases. 

We first assume that the plasma is hotter in the streaming direction, e.g. T y > T, that 
means the anisotropy is positive. Thus, the aperiodic solutions of Eq. (jHJ), which correspond 
to the cumulative filamentation-Weibel mode, are plotted with solid lines in Fig. [2] and Fig. 
[3] for two different values of the temperature anisotropies, A = 1 and A = 0.5, respectively. 
In Ref. {2] it was shown that the growth rates of the Weibel instability in a thermally 
anisotropic plasma with a bi-kappa distribution function decrease in comparison to those 
obtained for the case of a bi-Maxwellian electron distribution function. The Weibel growth 
rates are plotted here with dashed line. On the other hand, in counterstreaming plasmas 
with sufficiently large temperature anisotropies we observe the opposite effect: the growth 
rates of the cumulative filamentation-Weibel instability obtained for bi-kappa anisotropies 
increase and become much higher than those obtained for counterstreaming plasmas with 
bi-Maxwellian anisotropies. This is the case in Fig. [2] (a) and Fig. [3] (a), provided a 




FIG. 2: The growth rates of the filamentation instability (solid lines) obtained from ([8j) for two 
values of relative streaming velocity: (a) vq/c = 0.1 and (b) Vq/c = 0.05, and which cumulate the 
Weibel instability effect due to a bi- kappa temperature anisotropy, A = 1. The coordinates are 
scaled as W= uj; L /uj pe and K= kc/oj pe . For comparison, with dashed lines are shown the growth 
rates of simply Weibel modes (vq = 0). 




FIG. 3: The same as in Fig. O but for a smaller anisotropy A = 0.5. 

sufficiently large relative streaming velocity, v , is used. Such behaviour is typical for the 
simple filamentations instability, which develops in isotropic plasma counterstreams (A = 0). 
In this case, the filamentation instability growth rates are solutions of Eq. (TTU|) . and they are 
shown in Fig. H](a) and (b) for two different streaming velocities, vq/c = 0.1 and vq/c = 0.05, 
respectively. 

Otherwise, if the streaming velocity, v , is small the kappa growth rates are lower than 
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FIG. 4: The growth rates of the filamentation instability (A = 0) obtained from (|10p for two 
different values of relative streaming velocity: (a) vq/c = 0.1 and (b) vq/c = 0.05. On the ordinate 
and abscissa are W= uJi/ui pe and K= kc/u pe , respectively. 

the Maxwellian, and this is shown in Fig. ([2]) (b). 

Due to a finite temperature in both plasma counterstreams the purely growing solutions 
appear to resemble a Weibel regime and therefore their existence is limited to the small wave 
numbers less than a cutoff value, k < k c , which is the nontrivial solution of (jSJ) to the limit 
of uii = 



Finally, when the contrastreaming plasma is colder in the streaming direction, i.e., 
T y < T, the temperature anisotropy becomes negative, — 1 < A < 0. In this case, the 
effective velocity anisotropy imposed to the plasma particles by the relative motion of the 
counterstreaming plasmas, is reduced, and this explains the lower growth rates obtained 
in Fig. [5] for the negative temperature anisotropies. For a small streaming velocity, but 
a sufficiently large negative anisotropy, the instability is suppressed and instead natural 
damped modes occur. Further improvement of stabilization is observed for those modes 
corresponding to the larger kappa indices, which approach the Maxwellian modes. 

To summarize, we have investigated the properties of purely growing electromagnetic 
instabilities in an anisotropic counterstreaming plasma and propagates perpendicular to the 
streaming direction. Specifically, we have considered for the first time the filamentation 
instability in counterstreaming plasmas with isotropic kappa distributions, and derived the 
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FIG. 5: The growth rates of the filamentation instability (solid lines) obtained from (JSj) for negative 
temperature anisotropies: (a) vq/c = 0.1, A = —0.5, (b) vq/c = 0.1, A = —0.95, and (c) vq/c = 
0.07, A = —0.5. On the ordinate and abscissa are W= LOi/co pe and K= kc/u> pe , respectively. 

dispersion relation ( jTUl) . and displayed its numerical solutions in Fig. [51 We found that 
the aperiodic instability can reach maximum growth rates significantly larger than those 
calculated for the plasma with a Maxwellian distribution function. This is contrary to the 
behaviour of the Weibel instability in a bi-kappa plasma, which has growth rates smaller 
than those obtained for a bi-Maxwellian plasma. 

Moreover, we have generalized the theoretical approach assuming that each of the coun- 
terstreaming plasmas exhibits a temperature anisotropy of a non-Maxwellian type, namely, a 
bi-kappa distribution function. In this case, the instability is driven equally by the free energy 
stored in the relative motion of counterstreaming plasmas and the temperature anisotropy 
of the plasma particles. Comparing Figures [2] and [3] with Figure |5j we observe that the 
filamentation instability growth rate is enhanced in the presence of the positive tempera- 
ture anisotropies, when plasma is hotter in the streaming direction. Higher growth rates 
are obtained for larger temperature anisotropies. This is the first effect of the so-called 
cumulative filamentation- Weibel instability, and it improves significantly the efficiency of 
the magnetic field generation, and provides further support for the potential role of the 
Weibel-type instabilities in the fast magnetization scenarios, e.g. gamma-ray bursts, etc. 

For negative temperature anisotropies the effective velocity anisotropy is reduced, and 
therefore the growth rates of the filamentation instability decrease, but they remain larger 
than those for the bi-Maxwellian plasmas. The unstable modes with largest kappa indices 
can be evenly suppressed. This is the second cumulative effect of the counterstreaming 
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plasmas and the (negative) temperature anisotropy, which contributes in this case to the 
stabilization of the plasma system. 
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